Fermionic and bosonic ac conductivities at strong disorder 
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We study the ac conduction in a system of fermions or bosons strongly localised in a disordered 
array of sites with short-range interactions at frequencies larger than the intersite tunnelling but 
smaller than the characteristic fluctuation of the on-site energy. While the main contribution cyo{u)) 
to the conductivity comes from local dipole-type excitations on close pairs of sites, coherent processes 
on three or more sites lead to an interference correction 0-1(0;), which depends on the statistics of the 
charge carriers and can be suppressed by magnetic field. For bosons the correction is always positive, 
while for fermions it can be positive or negative depending on whether the conduction is dominated 
by effective single-particle or single-hole processes. We calculate the conductivity explicitly assuming 
a constant density of states of single-site excitations. Independently of the statistics, cro(ij;) = const. 
For bosons cri(a;) oc log(C/a;). For fermions (Ti(cj) oc log[max(y4,a;)/a;] — log[max(_B,cj)/a;], where 
the first and the second term are respectively the particle and hole contributions, A and B being 
the particle and hole energy cutoffs. The ac magnetoresistance has the same sign as a\{u}). 

PACS numbers: 71.10.-w, 77.22. Gm, 74.25. N, 67.85.-d 



Interestingly, interference and many-body phenomena 
in strongly disordered insulators may be comparable to 
or even greater than those in clean materials. Long ago 
fermionic many-body cotunnelling at strong disorder has 
been addressed in the context of variable-range hopping 
in semiconductors The advent of modern cold-atom 
techniques and recent experiments on disordered su- 
perconductive films 0-l5j offered ample new opportuni- 
ties to observe interference and many-body effects at 
strong disorder, yet access them in both fermionic and 
bosonic systems. Recently it was shown in Ref. and 
later discussed in Q that interference effects in bosonic 
systems are particularly strong due to the constructive in- 
terference between all low-energy single- and many-boson 
cotunnelling processes, which leads, e.g., to a huge posi- 
tive magnetoresistance [1, Q and broadens the supcrfluid 
phase of the bosons. 

dc Transport at strong disorder requires inelastic hop- 
ping of charge carriers through the whole sample, as- 
sisted by absorption/emission of neutral excitations, e.g., 
phononsfiol. \\V^. ac Conductivity may come from local 
excitations and does not vanish even in the limit of in- 
finitesimal dissipation [12[. The simplest local excitation 
is a particle-antiparticlc dipole created by external field 
on a pair of neighbouring sites [l3|- impurities or quasi- 
localised states. To study ac ma gne toresistance and the 
Hall effect in hopping insulators|14l 4l'/J one has to con- 
sider processes on three or more sites|13|. So far stud- 
ies of magnetotransport and interference phenomena con- 
cerned fermionic charge carriers (electrons) and did not 
address the influence of quantum statistics. 

In this paper we calculate the ac conductivity of 
strongly disordered fermionic and bosonic insulators at 



zero temperature. We demonstrate that coherent pro- 
cesses on three or more sites cause the interference cor- 
rection to the conductivity, that depends on the statis- 
tics of the charge carriers. For bosons the correction is 
always positive, while for fermions it may be positive or 
negative if the conductivity is dominated by particles or 
holes respectively. The latter holds even in the case of 
non-interacting fermions. We calculate the conductivity 
explicitly for sufficiently large frequencies assuming that 
the density of the charge carrier states is constant at re- 
spective energies. 

Our results apply equally to charged particles in disor- 
dered media and to neutral cold atoms in optical lattices. 
In the latter case the conductivity should be understood 
as a response to a tilt of the external potential. 

Model. Strongly disordered insulator may be modelled 
as an array of random-energy sites, e.g., electron puddles, 
grains, impurities, quasi-localised states, with a weak in- 
tersite tunnelling. We assume that the tunnelling con- 
serves particle spin and that the characteristic on-site in- 
teractions significantly exceed the characteristic energies 
of conducting excitations. Then one may consider spin- 
less particles, multiplying the conductivity by the spin 
degeneracy in the end. The generic Hamiltonian of the 
charge carriers reads 



(1) 



where 6j and br are the (bosonic or fermionic) parti- 
cle creation and annihilation operators on site r, trr'- 
the intersite tunnelling clement; Ur{nr) is the energy of 
rir = 6j6r particles on site r- a large random function of 
rir, which accounts for the interaction of particles and a 



2 



random potential. In particular, these energies prevent 
superfluidity in the case of bosonic particles. The dis- 
order manifests itself in the fluctuations of the on-site 
energies, Ur{nr), intersite couplings irr', and random lo- 
cations of the sites. 

The intersite interactions are assumed to be negligible, 
which corresponds to a system of neutral cold atoms a 
Josephson network with large self-capacitances [Iq . and 
can be also justified in systems with weak Coulomb in- 
teractions, provided all relevant energies sufficiently ex- 
ceed the Coulomb gap[ll|. We assume that each occu- 
pation number rir uniquely determines the many-particle 
state on site r because of sufficiently strong on-site inter- 
actions, leading to large energy gaps between different 
orbital states of the particles. Eq. ([T} is a generic 
Hamiltonian that describes, for instance, the Anderson 
model or the disordered Bose-Hubbard model in the case 
of fermions or bosons respectively. 

The strongly insulating regime requires the tunnelling 
trr' to be small compared to the fluctuations of the on- 
site energies and to sufficiently quickly decay with dis- 
tance. Thus, the groundstatc of the Hamiltonian ([T]) is 
close to the groundstatc of the on-site energies, that cor- 
responds to the conflguration of integers minimising 
the flrst term in the right-hand side of Eq. ([I]) @. At 
low temperatures and frequencies only the lowest-energy 
excitations are important on each siteQ, corresponding 
to the occupation numbers -I- 1 and — 1, which will 
be referred to as particle and hole excitations on site r 
respectively, both in the case of bosons and fermions. 

In the case of fermions the important lowest-energy 
excitations correspond to only one orbital state on each 
site, the interaction with the other electrons being equiv- 
alent to a static potential, and, without loss of generality, 
we may assume that the occupation numbers take only 
values ~ and — 1. In the case of bosons we as- 
sume ^ 1 for simplicity; however, the conclusions of 
this paper hold at arbitrary average characteristic occu- 
pation numbers. 

We calculate the conductivity assuming that the exci- 
tation density of states v per site is constant up to the 
characteristic cutoffs and E'^ of the particle and hole 
excitation energies, and that the frequency exceeds the 
characteristic matrix elements of the intersite couplings 
but is smaller than the characteristic excitation energies, 

Jrr' < < iy-\max{EP,E''), (2) 

where Jrr' = ^rr' for fermionic charge carriers, while for 

bosons Jrr' = ^rr' ('^r^'r' ) ■ 

A uniform external field E{lu) induces a current, 
Ii2{uj), on each pair of sites 1 and 2, which at low temper- 
atures can be expressed through the retarded correlator 
of Ii2 and currents Irr' on all the other pairs of sites using 
the Kubo formula: 



n 



Il2{uj) ^ i{2uj) ^y^ni2.rr' {^^)Ur: 



(3) 
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Here /rr' = iq [t.r'rbl,br — trr'blbr'j and Urr' = (r' — r)E 
are the current operator and the effective voltage drop 
between sites r ad r', q- the particle charge, E- the uni- 
form amplitude of the electric field. In absence of exter- 
nal magnetic field the average current density evaluates 



2 



lEl 



r), 



(5) 



n being the concentration of the sites, (. . .)dis is our con- 
vention for the disorder averaging. 

Two-site conduction. At infinitesimal dissipation 
the conduction comes from the resonant absorption of 
the electromagnetic field quanta uj by local particle- 
antiparticle excitations on sparse pairs of sites |13|. Cre- 
ating more complicated excitations is suppressed by the 
smallness of the tunnelling. 



a 



{Jrr')d^sdr' ^1. (6) 



Clearly, the particle hopping between two sites is not af- 
fected by quantum statistics, except maybe for the value 
of the coefficient before the hopping rate. We evaluate 
first the contribution of these trivial two-site processes 
to the conductivity before analysing the interference and 
many-body corrections to it. We consider the limit of 
low temperatures T <^ uj. 

Two-site processes correspond to the terms with r — 
1(2), r' = 2(1) in Eq. dS]). Evaluating the correlator of 
the current 1 12 with itself we obtain the real part of the 
two-site conductance: 



Gl2=TTq^0J-^\Jl2\^S{E2 



El-Cj), 



(7) 



where Ei and E2 arc the energies of the particle on the re- 
spective sites, the ground state being close to the particle 
residing on site 1. The delta- function in Eq. ([7]) rcfiects 
the resonant absorption. Actually, in presence of a finite 
weak dissipation it should be ascribed a certain width 7, 
the concentration of resonant pairs of sites ~ nji/ being 
very small. Unless quenched disorder is long-correlated, 
the intersite couplings J12 = J|r2-ri| and the on-site en- 
ergies £'1^2 fluctuate independently. 

The conductances ^ between pairs of sites lead to the 
conductivity 



2 2 

(To ~ Trn q 



(8) 



Vdip{^,0 = 2((S(i?2 - El - uj))dis being the density of 
states of dipole excitations, particle-hole pairs of size ^, 
factor 2 accounting for the 2 possible polarisations of a 
dipole at the same Ei and E2. 
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For a constant density of states v of the single-site 
particle and hole excitations, considered in this paper, we 
find v^iy = 2i/^w, which yields a frequency-independent 
conductivity 
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(9) 



Let us emphasise that the conductivity tTo(w) is con- 
stant only in the frequency range under consideration. 
At higher frequencies, violating the second of inequali- 
ties ([2]), the conductivity decays due to the lack of suffi- 
ciently high-energy states. At smaller frequencies, when 
the first of inequalities ([2]) no longer holds, the conduc- 
tivity is described by the famous Mott's formula [l^ [l9j 
and decreases (x uP' . 

Frequency-independent conductivity, Eq. ([9]), can be 
understood from the Mott's formula as follows. A small 
voltage F, applied to a resonant pair of sites, Ei — E2 = 
w, with a small intersite coupling J <^ uj, makes a pertur- 
bation with the off-diagonal entry ^ qVJ/{Ei — E2) oc 
uj~^. The Mott's formula, on the opposite, applies when 
the typical coupling exceeds the frequency, uj < {J)dis- 
Then the matrix element of the intersite transitions is 
non-perturbative in J/uj and does not depend on fre- 
quency (except, maybe, for a logarithmic factor). Be- 
cause the conductivity is quadratic in the matrix element, 
ct(w) oc aMott{'jj)<^~'^ = const at high frequencies under 
consideration, uj :$> {J)dis- 

Three close sites is the smallest cluster of sites that 
accounts for the interference and non-trivial quantum 
statistics effects. The contribution of larger clusters to 
the conductivity is suppressed due to the smallness of the 
intersite tunnelling, Eq. ([6]). 

There are two possibilities for the lowest-energy single- 
site excitations in a cluster of three sites: 1) a hole can 
occur on one site and an extra particle- on each of the 
other two, or, vice versa, 2) a particle excitation can oc- 
cur on one site and a hole- on each of the other two. 

In the first case the dynamics of the three sites is equiv- 
alent to the hopping of a single particle between these 
sites, described, both in the case of bosons and fcrmions, 
by the effective Hamiltonian 
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where E^ is the energy of the particle on site r, Ei < 
E2,E3. The three-site Hamiltonian in the form (ITOl) has 

II) to 



been used in a number of works (cf. e.g. [3, IE 
study the Hall effect in hopping insulators. 

Single-hole hopping. In the second case the dynamics 
of the three sites is equivalent to the hopping of a lack 
of a particle between these sites. However, the effective 
Hamiltonian will depend dramatically on the statistics of 
the charge carriers. 



Indeed, the state of several fermions on several sites is 
antisymmetric with respect to the permutations of the 
fermions, which makes the sign of the tunnelling element 
between two sites depend on the occupation numbers of 
the other sites, {hOk\blbk\0^lk) = (-l)^^=='+i for i < 
k — 2 [2l|. For bosons all the signs are the same. From 
Eq. ([T]) we find the effective Hamiltonian which describes 
the hopping of a hole on three sites: 



TjhoLe 
^3 



El — J21 TJ31 
— J12 E2 —J32 
T'h3 —J23 E3 



(11) 



where the upper and lower signs apply to bosons and 
fermions respectively, Er is the energy of a state with a 
hole on site r, Ei < E2,E^. 

The hopping of a single particle or hole is not a many- 
body problem effectively. Nevertheless, Eqs. PH)) and 
PT|) show that the parameters of the respective single- 
body Hamitonian qualitatively depend on the statistics: 
for bosons (particles and holes) and fermionic particles 
the signs of all tunnelling elements are the same, while 
for a fermionic hole the sign is alternating. Below we 
demonstrate that this difference manifests itself, under 
certain conditions, in the sign of the interference correc- 
tion to the ac conductivity. 

The currents between three sites may be found, using 
Eqs. ([3]) and (jl]), straightforwardly from the Hamiltoni- 
ans (fTU]) and pT|) . Because the interference correction to 
the conductivity has a relative smallness a, Eq. (I6|), com- 
pared to the two-site contribution ([9]), it is sufficient to 
find the currents up to the third order in the small cou- 
plings J, corresponding to the first-order-in-a correction 
to the two-site current I12 = G12U12, Eq. ([7]). 

Again, conduction requires some resonant excitation 
be present on the three sites. Below we assume that sites 
1 and 2 in the cluster under consideration are resonant, 
i.e. £'2 — -Bi ~ w. The probability of finding one more 
resonant pair in the same cluster ^ is negligible. 

In absence of external magnetic field all the couplings 
may be chosen real and positive, Jrr' = \Jrr'\- Diagonal- 
ising the Hamiltonians (fTO|) and (fTT]) , and evaluating the 
correlators of the currents, from Eqs. ([3]) and (|4]) we find 
the currents due to the hopping of particles and holes on 
three sites: 



h2^Trq'L0~'S{E2-Ei-Uj)-> \J12VU12 



± J12 J23>/l3 



C/13 



32 



Ez — E2 E^ — El 



Il3 = ±'KqUJ S{E2 ~ El - Uj)Ji2J23Jl3 

I23 = ±7r(7^a;~^(5(i?2 ~ Ei ~ w) J12 J23-^i3 



E3 - E2 

U12 
E3 — El 



(12) 

,(13) 
■,(14) 
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where the upper signs apply to fermionic particles and to 
bosonic particles and holes, the lower signs- to fermionic 
holes only. 

The first line of Eq. ([T^ is the current due to the 
processes on the two resonant sites, while the second line 
of Eq. (in]) and Eqs. (US]) and (HI]) are respectively the 
interference correction to the two-site current and the 
extra currents due to the presence of the third site. 

Bosonic interference correction. In most three-site 
clusters under consideration, the excitation energy 
on the third site greatly exceeds the energies Ei,E2 ^ oj 
on the two resonant sites due to the condition ([2]) . 

In the case of bosons the direct current between sites 
1 and 2 is increased by the presence of the third site at 
-©3 3> El, E2, yet the extra current /13 ss 1^2 7 flowing 
through the third site, additionally enhances transport 
between sites 1 and 2. This reflects the general principle 
that all low-energy cotunnelling processes in a bosonic 
system interfere constructively, effectively enhancing low- 
energy transport 0, HI] . 

To find the conductivity it is convenient to average 
first the currents in a three-site cluster with respect to 
the energy of the excitation on site 3. The positions of 
sites fiuctuate independently of the couplings and can be 
averaged separately. Using Eqs. (O, (fT^ - (fH)) wc find the 
interference correction to the conductivity ctq, Eq. ([5]), 
for a bosonic system with a constant density of states of 
single-site excitations at w ^ EP-'^: 

cri (cj) = A log{EP /oj) + A log{E''/oj) (15) 

A = 2iTd-^n^y^ j (JcJeJ|g_ci)*.^'c?V'C- (16) 

The logarithms in Eq. ()15p come from integrating 
Eqs. (in])- (Hi]) wrt the third-site energies at w <C £'3 <C 
E^^^ and must be smaller than to ensure the valid- 
ity of the perturbation theory. The frequency lo serves 
as an effective low-energy cutoff due to the saturation of 
the average current corrections at small energies on the 
third site, Q < E^ — Ei < E2 - Ei = uj. If the frequency 
w exceeds the energy cutoffs E^ or of the particle or 
hole excitations, the respective contribution in Eq. (ITS)) 
vanishes. 

Eqs. (IT5)) and show that in a bosonic system both 
particle and hole contributions to the interference correc- 
tion to the conductivity are positive. 

Fermionic interference correction. The single-particle 
and single-hole processes contribute to the interference 
correction with different signs in the case of fermions: 

fTi = yllog[max(£;P,a;)/w] - Alog[max(£''', w)/a;], (17) 

where constant A is defined by Eq. (ir6| . The cutoffs E^ 
and E^ depend on the level of doping of the insulator. 
If the band of the localised states is nearly empty, then 
EP ^ E^, and the interference correction to the conduc- 
tivity is positive. An almost filled band corresponds to 
EP ^ E^, leading to a negative correction (|17p . 



In the whole frequency range. Eq. the correc- 
tion is positive or negative if the conduction is domi- 
nated respectively by particles or holes. At sufficiently 
small frequencies, w <^ E^'^, the correction is frequency- 
independent, (Ti = A\og{EP /E'^). Despite this, it can 
be separated from the trivial two-site contribution (Tq 
using the suppression of interference by magnetic field, 
as we discuss below. If the frequency lies between the 
two cutoffs, EP''^ -C a; <C E'^-p, the conduction is dom- 
inated by the excitations with the higher cutoff, cti = 
±A\og{EP''^ /oj). Let us notice, that in the latter regime 
of an effectively one type of fermionic charge carriers, the 
frequency dependency of ai is the same as for bosons, 
Eq. (IT5]) . but has a different sign if the charge carriers 
are holes. 

Suppression by gauge fields. In presence of magnetic 
field the tunnelling elements of charged particles acquire 
phases 

Jrr' = I Jrr' | CXp (^tqcr' ^ A(^)d^ ^ , (18) 

where A is the vector potential, related to the magnetic 
field, and we neglected the modification of the on-site 
wavefunctions. In the case of neutral atoms in optical 
lattices, the role of magnetic field may be played by artifi- 
cial gauge fields ^induced by rotation or the Berry phases 
of atomic levels [2|. 

The phase factors fluctuate randomly, due to the ran- 
dom relative positions of the sites. This does not affect 
the two-site contribution to the conductivity ctq but de- 
stroys the interference corrections to it if the magnetic 
flux through a characteristic three-site cluster exceeds 
one flux quantum qBX^ 3> 1, where A is the charac- 
teristic radius of the coupling Jrr' ■ 

Charge carriers with spin projections Sz acquire addi- 
tional energies —s^B in magnetic field, which has a neg- 
ligible effect on the conductivity provided \szB\ ^ EP^'''' . 
Thus, the interference correction determines the sign of 
the magnetoresistance: positive for bosons, and positive 
or negative for fermions depending on the doping level. 

Discussion. Wc studied ac transport in strongly dis- 
ordered systems. The conductivity is dominated by 
statistics-independent processes on pairs of sites, while 
larger site clusters give rise to the interference correc- 
tion, which is always positive for bosons and can have 
either sign for fermions. 

Indeed, different commutation rules for bosonic and 
fermionic variables lead to different character of interfer- 
ence effects [2^. All low-energy bosonic processes inter- 
fere constructively, effectively enhancing transport Q and 
leading to a positive bosonic interference correction. 

Fermionic many-body processes can have different 
signs. We have shown that particle-dominated fermionic 
processes give a positive correction, hole-dominatcd- 
negative, even in the case of non-interacting fermions. 
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At first glance, transport of non-interacting particles and 
holes is a single-body problem and should not depend on 
the statistics of the charge carriers. Nevertheless, the 
parameters of the effective single-body Hamiltonians de- 
pend on the statistics, which in the case of fermions mani- 
fests itself in the alternating signs of the effective intersite 
couplings, Eq. pT|) . 

We have calculated the conductivity explicitly, 
Eqs. (|9]), (fTS]) and ([T7| . under the assumptions of a con- 
stant density of states, low temperature, and negligible 
intersite interactions at relevant energies. However, our 
qualitative conclusions for the sign of the ac magnetore- 
sistance are valid for arbitrary interactions, densities of 
states, and temperatures smaller than the on-site energy 
fluctuations. Experimentally one can suppress the in- 
terference correction by magnetic field or artificial gauge 
fields and thus verify its sign. In a strong fermionic insu- 
lator the interference contribution to the ac magnetore- 
sistance changes sign when changing the doping level. 

Our results can be tested straightforwardly in 
experiments on spin-polarised fermionic [23| and 
bosonic(23| cold atoms localised in incommensurate 
optical lattices 25 1 or in a random potential [1^ H^l- 
The interactions in these systems are short-ranged, 
yet the disorder and the interaction strength are easily 
controlled 0. These systems may be in the insulating 
state, which implies the smallness of the intersite 
coupling compared to the fiuctuations of the on-site 
energies. The respective frequencies, cf. Eq. ([2]), lie in an 
easily accessible kilohertz range and exceed the typical 
temperatures T < lOnK. Thus, all the assumptions 
of this paper are fulfilled in such experiments, which 
allows one to verify not only our qualitative conclusions 
for the signs of magnetoresistance but also the explicit 
dependencies of the conductivity on frequency. 

The discussed effects can be observed also in disordered 
superconductive films in the insulating state. Again, 
strong insulation indicates of the smallness of the tun- 
nelling. The characteristic frequencies used in the re- 
cent experiments, Refs. (2? " 



23, w ~ 10 



lOOGHz, sig- 
nificantly exceed achievable temperatures T .IK ^ 
IGHz. The behaviour of the ac conductivity in mag- 
netic field would help one to identify which charge car- 
riers, electrons or Cooper pairs, dominate transport in 
those materials. 
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